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ABSTRACT 

There is increasing evidence that some massive globular clusters (GCs) host multiple 
stellar populations having different heavy element abundances enriched by supernovae. 
They usually accompany multiple red giant branches (RGBs) in the color-magnitude 
diagrams (CMDs), and are distinguished from most of the other GCs which display 
(-h \ variations only in light element abundances. In order to investigate the star formation 

Q^f histories of these peculiar GCs, we have constructed synthetic CMDs for u Cen, M22, 

and NGC 1851. Our models are based on the updated versions of Yonsei- Yale (Y 2 ) 
isochrones and horizontal branch (HB) evolutionary tracks which include the cases of 
enhancements in both helium and the total CNO abundances. To estimate ages and 
helium abundances of subpopulations in each GC, we have compared our models with 
the observations on the Hess diagram by employing a \ 2 minimization technique. We 
find that metal-rich subpopulations in each of these GCs are also enhanced in helium 
abundance, and the age differences between the metal-rich and metal-poor subpopula- 
tions are fairly small (~0.3— 1.7 Gyr), even in the models with the observed variations 
in the total CNO content. These are required to simultaneously reproduce the observed 
extended HB and the splits on the main sequence, subgiant branch, and RGB. Our 
results are consistent with the hypothesis that these GCs are the relics of more massive 
primeval dwarf galaxies that merged and disrupted to form the proto-Galaxy. 

Subject headings: globular clusters: individual (w Centauri, M22, NGC 1851) — Galaxy: 
formation — stars: abundances — stars: evolution — stars: horizontal-branch 



1. INTRODUCTION 

During the past decade, new observations have drastically changed our conventional view 
on the stellar population content of globular clusters (GCs) in the Milky Way. It is now generally 
believed that most GCs possess two or more stellar populations rather than a single population hav- 
ing the same age and homogeneous chemical composition. Photometric observations have revealed 
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magnitude diagram (CMD; e.g.. Lee et al. 


„^ „ „ ^ — 
1999: Pancino et al. 2000: Bedin et al. 2004: Sieeel et al. 


2007: Piotto et al. 20071: Milone et al. 


200 


5: 


Marino et al. 


2008: 


Piotto 


200< 


3: Moretti et al. 2009: 


J.-W. Lee et al. 


2009a: 


Ferraro et alJbood 


Han et al. 


2009 


: Roh et al. 


201V 


. Spectroscopic stud- 



ies have establis hed star-to-star variat ions in light element abundances s uch as C, N, O, Na, Mg, 
and/or Al (e.g., Carretta et al. 2009al jbl. and references therein; see also Gratton et al. 2012a . for 



D'Antona k Caloi 


2004: 


Norris 


20041: 1 


Lee et al. 2005. 2007: Piotto et al. 2005. 


2007; 


Renzini 2008: 


Yoon et al. 2008: 


DuDree et al.1 


2011: 


iMoni Bidin et al.l 


2011: 


Moehler et al. 


201 J: 


Gratton et al. 


2012bh. As proposed bv 


Carretta et al. ( 


2010bh and 


Marino et al. 


( 


201 lb!) . now "normal" GCs can 



be defined, in terms of chemical composition, as those that have homogeneous abundance in heavy 
elements, such as calcium and iron, but show variations in light elements such as a Na-0 anticor- 
relation. The presence of multiple populations and chemical inhomogeneity in thes e GCs is mainly 
consi dered to be due to the pollution from the fast-rotating massive stars (FRMSs 



Decressin et al 



2007) and /or intermediate mass asymptotic giant branch (IMAGB) stars (jD'Antona &: Ventura 
20071 : 1 Ventura fc D'AntonJl2008h . 



However, some massive GCs, such 


as oj Centauri (Lee et al. 1999: Bcdin et al. 2004) 


, M54 


(Lavden & Saraiedini 


1997: Sieeel et al. 


20071). M22 (Marino et al. 


2009: Da Costa et al. 


2009: 


J.-W. Lee et al. 2009a 


), NGC 1851 


Han et al. 2009: Carretta et al. 2010a), Terzan 5 (Ferraro et al. 


2009), and NGC 2419 


(Cohen et al. 


2010; 


di Criscienzo et al. 2011). are distinguished from the nor- 



mal GCs. They show evidence of supernovae (SNe) enrichment in addition to the spreads in light 
elements. The discrete distributions of red giant branches (RGBs) observed in the CMDs of these 
GCs, together with spectroscopic variations in heavy element abundances, indicate successive for- 
mation of stellar generations from the gas enriched by SNe. They are generally thought to be 
the relics of mo r e massive primeval dwarf galaxies that disrupted by and merged with the G alaxy 
(|Lee et al.lll999l : lGnedin et al.ll2002l : iBekki fc Freemanll20oi : iRomano et al.ll2007l : iBokerlbood ). and 
therefore, have important implications on the hierarchical merging paradigm of galaxy formation. 

The purpose of this paper is to investigate the star formation histories of three GCs that belong 
to the latter category, i.e., u Cen, M22, and NGC 1851. In order to estimate ages, metallicities, and 
helium abundances of subpopulations in these GCs, we have constructed synthetic CMDs which 
best reproduce the observations. The construction of synthetic CMDs were performed in a self- 
consistent manner from main sequence (MS) to horizontal branch (HB), whose morphologies vary 
with age, metallicity, and helium abundance. In the following section, we describe our construction 
of stellar population models. The synthetic CMDs, together with our best estimates of model 
parameters, for the three GCs are given in Section 3. Finally, based on our results, we briefly 
discuss possible formation scenarios of these GCs in Section 4. 
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MODEL CONSTRUCTION 



2.1. Stellar Population Models 



also 



The synthetic CMDs were constructed following t he technique s deve loped bv lLee et al.l (|l990l . 



1994) for the HB and RR Lyrae varia ble stars, and bv lPark k Led (|1997l ) for the MS to RGB (see 



Lee et al.ll2005l ; lYoon et al.l 12008). Our models are based on the latest version of Yonsei-Yale 



(Y 2 ) isochrones and a ne w set of Y 2 HB evolutionary tracks for the normal and enhanced helium 



abundances, respectively (|Yi et al.l 1200a : Y.-W. Lee et al. 2012, in preparation) . The HB evolu 



tionary tra cks used in our modeling; are the updated version of those presented by lLee &: Demarque 



(jl990l ) and lYi et al.1 (119971 ). and they have been calculated using the same code and input physics 



employed in the construction of Y 2 isochrones. Recent observations have shown that [CNO/Fe] 



Yong et al. 


2009; 


Marino et al. 


2011b. 


2012alJbl: 1 



order to reflect this in our modeling, we have expanded the parameter space of Y 2 isochrones and 
HB evolutionary tracks to include the cases of normal and enhanced nitrogen abundances ([N/Fe] 
= 0.0, 0.8, and 1.6; Y.-W. Lee et al. 2012, in preparation). The observed variations in the total 
CNO abundance were reproduced by interpolating these nitrogen enhanced isochrones and HB 
evolutionary tracks. Our test simulations with varying N and O abundances show that, once the 
total CNO sum ([CN O/Fe]) is held const a nt, both N and O have almost identical effect s on the 
HR diagram (see also lVentura et al.ll2009l : ISbordone et alJl201ll ; IVandenBerg et al.ll2012i ). Model 
ingredients adopted in our construction of synthetic CMDs are summarized in Table 1. 

In the evolution of low-mass stars, mass loss on the RGB plays an important ro le in determining 
the morphology of HB . In order to estimate the amount of this m ass-loss, followin g iLee et al.l (119941 ) 
and lRey et al.1 (120011 ). we have adopted the empirical formula of iReimerd (| 19771 ) and calibrated its 
coefficient, r/, to the HB morphologies of old GCs in the inner halo (Rgc < 8 kpc). We have first 
integrated mass-loss rate with time steps along the evolutionary tracks from the low MS to the RGB 
tip, and tabulated the resulting values of mass-loss as functions of rj, age, and chemical composition. 
Using these tables, together with the values of the RGB tip masses (calculated without the mass- 
loss) taken from the Y 2 isochrones, we constructed a series of synthetic HB models to form HB 
isochrones. Figure 1 compares these HB isochrones with the inner halo GCs on the HB morphology- 
metallicity diagram, under different assumptions as to the value of rj. Assuming that these inne r 
halo GCs are coeval and have a mean age of 13 Gyr (jMarm-Franch et al.ll2009l ; iDotter et al.ll2010i ). 
we obtained r/ ~ 0.53, which best reproduces the observed correlation. We have used this same 
value of T) in all of our model constructions. For the populations enhanced in helium and CNO 
abundances, the same value of r\ was adopted under the assumption that the mass-loss efficiency 
parameter, rj, would not be changed with helium or CNO abundances. Note, however, that the 
changes in gravity and luminosity caused by helium and CNO enhancements were still reflected in 
our calculation of the amount of mass-loss. Note also that the empirical fitting factor, rj, depends 
on the stellar libraries chosen and the ages assumed. For example, rj would be increased to ~0.59, 
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if we adopted 12 Gyr, instead of 13 Gyr, as the mean age for the inner halo GCs. 



Theoretical HR diagrams (log L vs. log T e g) of our models were conver ted into obs e rvable 



Green et al.l 



Green et al 



1987) 



19871 ) 



quantities (magnitude vs. color) by empl oying the semiemp i rical c olor table of 
and the synthetic model atmospheres of ICastelli Kuruca (|2003l ). While the 
table has been adopte d to the Y 2 isochrones and pro ven to show reasonable agreements with ob- 
served CMDs of GCs (|Yi et al.ll200ll ; iKim et al.ll2002l ). it is onl y availab l e for t he Johnson-Cousins 
UBVRI system. Therefore, foll owing the method describ ed in lYi et al .1 (|1995l ). we have generated 
synthetic color ta bles from the ICastelli fe: Kuruca (|2003l ) stellar spectral library for the Johnson- 



Cousins UBV RI (lBesselllll990ah. Hubble Space Telescope (HST) AC S/WFC, HST WFC3/UVIS, 



and Ca-uvby ([Crawford fc Barneslll970l : lAnthonv-Twarog et al.lll991f) fil ter systems. The bolomet 



rop oints of each filter sy stem (i.e.. iBessellll 1 990bl for the Johnson-Cousins UBVRI, 



ric corrections were computed using the definition of lBessell et al .1 (119981). and t he magnitudes and 
color indices were calibrated with the model flux of Vega (|Castelli fc Kurucahoosh to the relevant 



zc- 



Bedin et al 



2005 



and Potter et al.l 120071 for the HST ACS/WFC and HST WFC3/UVIS, and lCrawford fc Barnes 
1970 for the Stromgren uvby filter systems). The hk index was separately normalized with the 
model flux of the standard star HP 83373, which has the closest spectral type to that of Vega in 

4.0, [M/H] = 0.0, 



the catalog of I Anthony- Twarog et al 



and fturb = 2.0 km/s (jSoon et al 



1993 



199lh. assuming T P g = 1 0550 K, log g 
Castelli & Kuruca bood ). 2 



Figure 2 compares our isochrones and zero-age HB (ZAHB) models transformed by using the 
two different color tables. They are for the chemical compositions and ages appropriate for the 
most metal-poor and the most metal-rich (also helium-rich) subpopulations in cj Cen (see §3.1 
below). While the models are in general agreement, there are minor but noticeable differences in 
the slope of the metal-rich RGB, especially in panel (d), i.e., B — R color. Assuming that these 
differences in the Johnson-Cousins UBVRI system would be similar in the adjacent passbands 
of the HST filter system, we have applied small corrections to the synthetic magnitudes of the 
HST ACS/WFC F435W and F625W passbands. For instance, for the F435W band, we added the 
differences in the B band magnitudes between the two tables, to the synthetic F435W magnitudes. 
Similarly, we used the differences in the R bands for the F625W bands. The corrected colors are 
compared with synthetic colors in panel (e). In case of the UV (F225W, F275W, and F336W) and 
the broad VI passbands (F606W and F814W) of the HST WFC3/UVIS and ACS/WFC systems, 
the differences between the two color transformations are relatively small (see Fig. 2a and 2c), and 
therefore, we did not apply the corrections to these passbands. The Stromgren b and y bands are 
very similar to the Johnson B and V bands, and because the two color transformations agree well 
for the metal-poor population in B — V color (see Fig. 2b), we also did not apply the corrections 
to these passbands. Inter stellar extinctions fo r the filter systems used were estimated according to 
Cardelli et al.l (119 89) and lSirianni et al.l (|2005l ). based on the E(B — V) values listed in the updated 



catalog of 



Harris 



4199J). 



"Adopted from Fiorella Castelli's web page, http://wwwuser.oat.ts.astro.it/castelli/colors.html 
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2.2. A Fitting Technique for Finding Best Parameters 

In order to fit our synthetic models to the observed CMDs and find the values of ages, metal- 
licities, and helium abundances of subpopulations in each GC, we have combined both eye fitting 
on the CMDs and x 2 minimization technique on the Hess diagrams. First, from the literature, we 
predetermined population ratio of subpopulations in a GC, metallicity ([Fe/H]) of each subpop illa- 
tion, and a difference in the total CNO content (A[CNO/Fe]) between subpopulations. Synthetic 
CMDs were then constructed for the relevant filter systems, and they were first compared with the 
observed CMDs by eye fitting. At this stage, initial guesses of age and helium abundance were 
assigned, and distance modulus, reddening, and [Fe/H] values of each subpopulation were fixed 
within the observational uncertainty. 

The best fitting values of age and helium abundance for all subpopulations were then si- 
multaneously determined by performing a x 2 minimization between the observed and synthetic 
Hess diagrams. This was in s pired by the techn i ques widely used i n dwarf galaxy studies (e.g., 



spi 



Harris & Zaritskyll200ll . 12004 : iTolstoy et al.ll2009l : iRubele et al.ll201ll ). For example, Figures 3 and 



4 illustrate our procedure for NGC 1851. Panel (a) in Figure 3 shows the CMD regions used in 
the x 2 minimization, and panels (b) and (d) display the Hess diagrams of the observed CMD. 
Using the predetermined parameters of each GC, we constructed a set of synthetic CMDs with an 
age resolution of 0.1 Gyr and helium abundance resolution of 0.01, including observational error 
simulations. 3 Based on these models, synthetic Hess diagrams were produced by averaging 1000 
simulations to minimize stochastic effects. The synthetic Hess diagrams were finally compared 
with the observed ones using the IDL 4 function XSQ_TEST, which performs the x 2 goodness-of-fit 
test between the observed and expected frequencies of a theoretical distribution. We have fixed 
helium abundance of the metal-poor first-generation (Gl) subpopulation to be normal value (i.e., 
Y = Y p + 2Z, Y p = 0.23), so that Gl was examined only for age variation, while the metal-rich 
second-generation (G2) subpopulation was investigated for the variations of both age and helium 
abundance, simultaneously. Panels (c) and (e) in Figure 3 are the model Hess diagrams with the 
minimum x 2 ■, an d Figure 4 is the distribution of reduced x 2 ix 2 ) values for the variations in age 
and helium abundance of each subpopulation. The 68% (1 a) confidence levels were calculated by 
bootstrapping with 1000 simulations for each subpopulation. The estimated 1 a errors are typically 
± 0.3 Gyr for the relative age and ± 0.02 for helium abundance (Y). Note that these errors are 
only for the case when the metallicities and the total CNO contents are held fixed, and therefore, 
the actual errors would be larger if the measurement errors of metallicity and [CNO/Fe] are not 
negligible. The uncertainties in the distance modulus and reddening, which are fixed in the analysis, 
as well as those in stellar models, are also not included in the quoted errors. Note, however, that 
they would have only limited impact on the estimation of relative differences in age and helium 
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abundance. 



3. SYNTHETIC COLOR MAGNITUDE DIAGRAMS 



3.1. oj Cen 



As numerous photometric and spectroscopic studies have been devoted to it, the most massive 
GC, to Cen, is indeed the most complex object among GCs in the Milky Way. It shows several promi- 
nent features in the CMDs such as a double MS, multiple subgiant branches (SGBs), multiple RGBs, 



elements (Lee et al. 


19991: 


Pancino et al. 200ol: Rev et al. 


2004 : Ferraro et al. 


2004: Bedin et al. 


2004: Sollima et al. 


2005al. 


200d. 20071 


Piotto et al. 2005: 1 


v^illanova et al. 20071: 


Bellini et al .1 


2010|; 


Freeman & Rodgcrs 


1975: 


Norris et al. 


1996: Suntzeff k KraftJll996l: Ijohnson & Pilachowskil 


2010; 


Carretta et al. 2010. 


a; Marino et al. 2011a). The multiple 


and discrete RGBs, together with their 



spectroscopic variations in the abundance of heavy elements, directly indicate that subpopulations 
in oj Cen were influenced by SNe. Meanwhile, the presence of a double MS and the fact that the 
bluer MS is more metal-rich than th e redder MS strongly suggest that there i s a large variation 



of he lium abundance in this cluster (IBedin et al .1 [2004J; iNorris! 120041 : iLee et al 



2005 



Piotto et al 



20051 ). The presence of unusually ex tended HB is also naturally explained by the subpopulations 
enhanced in helium (ILee et al.ll2005l ). 

Following the techniques and spirits employed in our previous investigation (|Lee et al. I l2005l ). wc 
have constructed new synthetic CMDs for to Cen. As described in §2.1, there have been two major 
improvements in our modeling. One is the adoption of the new Y 2 HB evolutionary tracks, together 
with the corresponding Y 2 isochrones, with which the enhancements in both helium abundance and 
[CNO/Fe] can be tested self-consistently. The other is the improvement of the color-temperature 
transformation including the "semi-empirical corrections" for the HST ACS/WFC F435W and 
F625W passbands (see Fig. 2e). Furthermore, the best fitting parameters are now determined with 
the method described in §2.2. In the construction of synthetic CMDs, the metallicities ([Fe/H]) 
for the five subpopulations were first adopted from the most r ecent spectroscopic observations of 



RGB stars (Johnson k Pilachowski 



earlier estimates (Norris et al. 1996; 



2010; 



Marino et al 



Suntzeff k Kraftlll996 



Rev et al.ll2004l : iSollima et al.ll2005al lbl iPiotto et al.ll2005l : IVillanova et al.ll2007 ; 



201 l aj), which are la r gely consistent with 



Smith et al 



2000 



Pancino et al. 


2002 


Johnson et al. 


2008 



2009). The dif ferences in [CNO/Fej among the five subpopulations were taken from the recent 
observations by iMarino et al.1 (I2012al ). and the adopted values are listed in table 2. We have then 



3 We adopted the Salpeter initial mass function, with the standard index (x = 1.35), in our construction of synthetic 
CMDs. The choice of x has only a negligible effect in our analysis, since we are dealing with the relatively bright 
parts of the CMD where the mass range is small in old stellar populations. 

4 IDL is the Interactive Data Language, a product of ITT Visual Information Solutions: |http:/ /www.ittvis.com 
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determined the best fitting ages and helium abundances by performing the x 2 minimization on 
the Hess diagrams employing the HST ACS/WFC F435W and F658N passbands. I n Figu res 5 
and 6, our synthetic models are compared with the observed CMDs by iBellini et al.l (|2010i ). and 
the best fitting parameters obtained by our model simulations are listed in Table 2. These model 
CMDs were specifically constructed to reproduce the observed CMDs, adopting the same HST 
filter systems, population ratios, and photometric errors. 

For the two most metal-poor subpopulations, which we refer to as Gl and G2, reasonable 
agreements with the observations, from MS to HB, are achieved with normal helium abundances 
(i.e., Y = Y p + 2Z, Y p = 0.23). The RR Lyrae stars belonging to these metal-poor subpopulations 
are also providing a very firm constraint (AY < 0.01) on the helium abundance as the luminosi- 
ties and periods of RR Lyraes depend very sensitively on helium abundance (see below). For the 
ages of these subpopulations, our fitting technique is suggesting 13.1 ± 0.2 and 13.0 ± 0.3 Gyrs, 
respectively, and therefore they could be coeval to within the uncertainty. For the intermediate 
metallicity subpopulation, G3, however, a large amount of helium enhancement (AY = 0.18 ± 
0.02) is suggested from our fitting technique to reproduce not only the bluer and fainter MS but 
also the EBHB stars. This is readily understood by the effects of helium abundance in stellar 
astrophysics. Because of the higher core temperature and decreased envelope opacity, helium-rich 
stars, in general, are bluer and brighter for a given mass. Since they evolve faster than helium-poor 
stars, helium-rich stars would have smaller masses for a given age, and therefore the helium-rich 
MS appears both bluer and fainter than the helium-poo r sequence on the isochrone (see, e.g., 



Cox &: Giulil Il968l : lib en <fe Faulknerl Il968l : iLee et al.1 120051 ) . For the same reason, helium-rich HB 



stars have small er total masses than do less helium-rich stars and end up on the blue HB or EBHB 



(|Lee et al.1 11994 ) . The best fitting age for this helium enhanced G3 is 12.0 ± 0.4 Gyr, which is 
about 1 Gyr younger than the reference subpopulations (Gl and G2). Similarly, for the two most 
metal-rich subpopulations, G4 and G5, our models require helium enhancements as much as that 
for G3, from the presence of EBHB stars and the round shape of the SGBs. If they were not 
enhanced in helium abundance, they should have red HB stars or red clumps and more flat tened 
SGBs at their metallicities, [Fe/H] > —1.0 (see Fig. 8 below, see also ISollima et al.ll2005bl ). For 



these subpopulations, our fitting technique is suggesting a common age of 11.4 ± 0.4 Gyr, and 
therefore, the age difference between Gl and G5 is predicted to be about 1.7 ± 0.5 Gyr. This 



a-elements and s-process elements (e.g., 


Norris & Da Costa 


1995; 


Busso et al. 


1999; 


Smith et al. 


2000; 


Johnson & Pilachowski 


2010; 


D'Antona et al. 


2011; 


Marino et alJhoila. 2012a: 


D'Orazi et al. 


2011; 


Valcarce & Catelan 


2011). Note that the variation in total CNO content between the sub- 



populations has the most significant impact on this age difference. As shown in Figure 7, if there 
were no variation in [CN O/Fe], our models suggest that all of the subpopulations in oj Cen wo uld 



be almost coeval (see also iMarm-Franch et al. 



2010 



D'Antona et al. 



2011 



Marino et al.ll2012al ). 



In order to illustrate, independently, the effects of metallicity, [CNO/Fe], age, and helium 
abundance on the overall morphology in CMD, Figure 8 compares the synthetic CMDs constructed 
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under different combinations of these parameters. The presented models are for the metallicities 
comparable to the most metal-poor and the most metal-rich subpopulations in oj Cen, Gl and G5. 
First, the AZ-only model in panel (a) is for the case in which G5 has enhanced metallicity but the 
same [CNO/Fe], age, and helium abundance with Gl. The AZ + ACNO + At model in panel (b) 
is for the case in which [CNO/Fe] is also enhanced and age is decreased for G5. Both enhance d 
[CNO/Fe] and decreased age are pushing the HB of G5 to even redder color (see lLee et al.lll994l ). 
while they have opposite effects on the MS turn-off (MSTO) and SGB regions (see also Fig. 7) 
so that no appreciable change is seen there. Finally, AZ + ACNO + At + AY model in panel 
(c) is for the case where the helium abundance of G5 is also significantly enhanced, as in our best 
simulation. As is clear from Figure 8, G5 would have the red HB in the models without helium 
enhancement, while the enhanced helium abundance in panel (c) pushes the HB stars into the 
EBHB region. It is also shown that the luminosity and color of MSTO, together with the SGB 
slope, are fairly affected by this variation in helium abundance. 

In our modeling, the effects of enhancements in helium and the total CNO abundances were 
only considered in the stellar interior, and it was assumed that their effects on stellar atmosphere 
would be negligible. For the effect of helium abun dance, we have confirmed this using the helium- 
enhanced model fluxes of lCastelli h Kuruca (|2003l ). 6 Figure 9a shows the model spectra of two stars 
with different helium abundances (Ay = 0.10) but the same temperature, gravity, and metallicity 
(i.e., T e g = 6000 K, log g = 4.5, and [Fe/H] = —1.5), which correspond roughly to stars near the 
MSTOs of G2 and G3 in u Cen. We can see that a large yariation in helium abundance (Ay = 0.10 ) 
has only negligible effect on stellar atmosphere (see also Girardi et al. 2007 : Sbordone et alJ[201ll ). 
In the UV region, for the wavelengths shorter than ~4000 A, however, some differences between 
the two synthetic spectra are noticed. Figure 9c shows that, in the worst situation, combinations 
of T e ff, log g, [Fe/H], and the helium variations of Ay = 0.10 can make differences in magnitude 
up to ~0.15 in the UV bands. In the case of the total CNO abu ndance, our assumptio n appears to 
be valid also in the optical passbands. Using synthetic spectra, ISbordone et al.l (|201ll ) have shown 
that the variations in C, N, O, and Na abundances mainly affect wavelengths shorter than ~4000 
A. They have also reported that the effect on the hk index is negligible (at most ~0.04 mag), while 
CMDs including U band could be affected by ~0.2-0.3 mag for the ~0.3 dex variation in [CNO/Fe]. 
Considering these uncertainties in the UV regime, the best fitting parameters in our models are 
determined based only on the optical bands and the hk index. 

Some HB stars in our models are in the R R Lyrae instabil ity strip, and they can provide 
additional constraints in the modeling. Following iLee et alj (119901), we have compute d the funda- 
mental periods (P/), using the equation of stellar pulsation (jvan Albada &: Bakerlll97ll ). for the RR 
Lyrae candidates in the synthetic HB. Figure 10 compares our synthetic stars to the observed RR 



The model atmospheres for enhanced helium abundances can be downloaded at Fiorella Castelli's web page, 
http://wwwuser.oat.ts.astro.it/castelli/grids.html They present the synthetic spectra for three metallicity grids. 
[Fe/H] = -1.5, -0.5, and 0.5, and two gravity grids, log g = 1.5 and 4.5, for AY = 0.10. For AY = 0.20, only 
metal-rich ([Fe/H] = 0.5) models are available. 
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Lyrae variables in oj Cen, where the m etallicities of 74 RR Lyrae stars are from the spectroscopic 
measurements by ISollima et al.l (|2006l ). In panel (a), the periods of c-type RR Lyrae stars are 
fundament al ized assuming the period ratio, P>,/P n ,h = 0.745, between the c-type and a6-type RR 
Lyare stars (jClement et al.ll200ll ; lAlcock et al.ll2000l ). In our models, only the two most metal-poor 
subpopulations, Gl and G2, have RR Lyrae stars. The HB stars from the three metal-rich and 
super-helium-rich subpopulations, G3, G4, and G5, are too hot to produce RR Lyare stars as they 
are on the EBHB region far beyond the blue edge of the instability strip. In panel (b), metallicities 
of the synthetic RR Lyrae stars have been assigned randomly according to a Gaussian distribution 
centered on the metallicities of the Gl and G2 (see Table 2), respectively, with (Type/H] = 0.13, 
which is roughly consistent with the observational error in [Fe/H] measurement. The predicted pe- 
riods, V magnitudes, and [Fe/H] from our model populations, Gl and G2, show good agreements 
with the observed RR Lyrae variables exce pt for the several sta rs having intermediate metallicity 
(—1.4 < [Fe/H] < —1.0). As pointed out bv lSollima et al.l (|2006l ). these stars are expect ed to have 



norm al helium abundance because their luminosities and periods are not increased (see iLee et al 



19941 ). This implies that a minor subpopulation having a similar metallicity but different helium 
abundance with G3 may coexist in u Cen. The presence of these normal- helium, intermediat e 
metallicity RR Lyrae stars, together with the HST WFC3 UV photometry by lBellini et al.l (|2010l ). 
suggests that oj Cen contains additional minor subpopulations between the five subpopulations 
adopted in this study. 



3.2. M22 



It has long been suspected that RGB stars in the massive GC M22 show a spread in heavy 
element abundance, but the presence of differential reddening, caused by its locat i on near the Galac- 
tic plane, has prevented a firm confirmation of its reality (e.g 



1982; 



Ivans et al 



Nor ris fe Freemanlll983l : iLehnert et al.lll991l : lAnthony-Twarog et al 



Hesser et al 



1977 



Pilachowski et al.l 



1995; Monaco et al. 2004 



20041 ). Recent observations, however, have shown that M22 indeed hosts two stellar sub- 



populations d i ffering in the abundance of heavy elements. Based on high resolution spectroscopy. 



Marino et al 



(2009, l2011bl ) identified two stellar groups separated by ~0.15 dex in [Fe/H], while 



Da Costa et al 



reported a somewhat larger difference in metallicity (~0.26 dex). Meanwhile, 



J.-W. Lee et al.l (|2009al ) has clearly shown that M22 has a discrete double RGB having different 
Ca abundance from the hk index of the Ca-by photometry. 

In the left panels of Figure 11, we show CMDs of M22 from S.-I. Han et al. (2012, in prepara- 
tion), which were observed with the Stromgren b, y and a new narrow band calcium filter specially 
designed to avoid a possible contamination from the adjacent CN bands. We c an see two well- 



separat ed RGBs in the hk index, which are very similar to those in Figure 1 of I J.-W. Lee et al 



(j2009al ). and a bimodal distribution of HB stars compose d of blue HB and EBHB. Based on these 
photometric data and the CMD from the ACS survey (jPiottol l2009T ) 5 . shown in Fi gure 12a, we 
constructed synthetic CMDs as described in §2. The difference in [CNO/Fe] between the two sub- 
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populati ons was adopted to be 0.13 dex from the recent observations (jMarino et al.ll2011bl . l2012bl . 
see also lAlves-Brito et al.ll2012l ). 



The right panels of Figures 11 and 12 are the final outcomes of our population models compared 
with the observed CMDs. The best fitting model parameters are listed in Table 3. The metal-poor 
subpopulation, which we refer to as the first-generation (Gl) population, contains the bluer RGB, 
brighter SGB, and blue HB stars, while the metal-rich subpopulation, which we refer to as the 
second-generation (G2) population, has the redder RGB, fainter SGB, and EBHB stars. In order 
to simultaneously reproduce these observed splits on the RGB and SGB, and the bimodal HB 
distribution, we need a small but important difference in metallicity (A[Fe/H] = 0.25 dex), and 
at the same time, a large difference in helium abundance (AY = 0.09 ± 0.04) between the two 
subp opulations. By appl ying a small difference in the total CNO content (A[CNO/Fe] = 0.13 
dex; iMarino et al.ll2012bl ). our models suggest that the best fitting age difference is very small 
(At = 0.3 ± 0.4 Gyr). The adopted dif ference in metallicity is similar to the observed difference 
obtained from Ca II triplet (~0.26 dex: lDa Costa et al.1120091). while it is somewhat larger than the 



difference obtained from [Fe/H] (~0.15 dex; 



Marino et al 



20091. hoilbh. This is probably due to a 



Marino et al. 



2009 



2011bl ) in the metal-rich 



small increase in calcium abundance (A[Ca/Fe] ~ 0.1; 
subpopulation, G2. O ur models are also consistent with the recent spectroscopic measurements by 
Marino et al.l (|2012bl ). which have shown that stars in the faint SGB are more metal-rich compared 
to the stars in the bright SGB. 

Figure 13 explains the effects of variations in the total CNO abundance, overall metallicity, and 
helium abundance, independently. If we assume only the difference in [CNO/ Fe] between the two 



subpopulations, to reproduce the split on the SGB, as suggested for NGC 1851 (jCassisi et al. 



2008 



Salaris et al.ll2008l ). our models fail to explain the presence of the RGB split in the hk index and the 
EBHB stars (the upper panels). Note that, for illustrative purpose, we have increased the difference 
to 0.5 dex from the observed value (~0.13 dex), because the observed difference is too small to 
reproduce the SGB split. Similarly, with the assumption of enhancement only in [Fe/H], our model 
for G2 also fails to reproduce the EBHB stars, while the split on the RGB can be explained (the 
middle panels). This is a natural consequence of the metallicity effect, as the HB morphology is 
getting redder with increasing metallicity (see, e.g. JLee et al.lll994l ). While this difference in [Fe/H] 
reproduces some split on the SGB, it is not sufficient to explain the observed difference. The bottom 
panels of Figure 13 show that additional enhancement in helium is needed to overcome the effect 
of metallicity on the HB and to reproduce the EBHB stars, as already explained in §3.1. Similarly 
to the case of u Cen, this enhancement of helium also helps to reproduce the observed split on the 
SGB (see Fig. 12d). 

The argument based on the population ratio further supports our hypothesis that the EBHB 
stars belong to the metal-rich subpopulation G2. The number ratio between the bluer and redder 



5 The photometric data are available at the web page of the ACS survey of Galactic GCs l|Saraiedini et al.ll2007 
Anderson et al. 20081 ) . http://www.astro.ufl.edu/~ata/public_hstgc/databases.html 
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RGBs in Fig u re 11 c and those in |j.-W. Lee et al.l (|2009ah is about 0.7:0.3. IPiottol (|2009|) and 



Marino et al.l (|2009l ) have also suggested a similar number ratio (0.62:0.38) between the bright 
and faint SGBs. Within the error, this is also consistent with the ratio of the EBHB stars which 
comprise roughly 25% of all HB stars in panel (a) of Figures 11 and 12. This would then suggest 
that the bright SGB, bluer RGB, and blue HB stars are all associated with the majority, metal-poor 
subpopulation (Gl), while the faint SGB and redder RGB stars are progenitors of the EBHB stars, 
associated with the minority metal-rich subpopulation (G2). 



3.3. NGC 1851 



Recent observations of NGC 1851, which is well-known for its unusual bimodal HB distribu- 



tion (e.g., ISaviane et al 



1998 



Walkerl Il998). indicate that it also belongs to the peculiar group 



of massive GCs influenced by SNe enrichment like u Cen and M22. Not only a double RGB has 



been revealed from several photometric observations (jJ.-W. Lee et al. 



Carretta et al. 



2009 



SB 



lHan et al. 



2009 



201 laT). but also some spreads in [Fe/H] have been reported from spectroscopic obser- 



vations (jCarretta et al.ll2010d . l2011bl ). NGC 1851 also show s other signatures fo r the presence of tw o 
subpopulations, such as a double SGB (jMilone et al.1 120081') . a CN bimodalitv (IHesser et al.1 11 982). 
and variations in ligh t and s-process element abundances (jYong &: Grunda hl 2008; Yong et al.ll2009l ; 
Villanova et al.ll2010i : ICarretta et al.ll2011bl ). The origin of the observed SGB split, together with 
its possible links to the RGB split and the bimodal HB distribution, however, is still under debate 
For instance, the SGB spil t can be explained by a diffe rence in age of ~ 1—1.5 Gyr (jMilone et al 



2008; Carretta et al. 2010c. 2011b 



by a factor of ^2-4 (ICassisi et al 



Gratton et al.l 12012 bT). by an increase of total CNO abundance 



20081 : ISalaris et al.l 12008 : lYong et al J 12009 : IVentura et all 12003 : 



D'Antona et al.ll20 09). or by a combined effec t of the differenc es in overall metallicity (A [Fe/H] 
0.15 dex) and helium abundance (AY 0.05; lHan et al.ll2009i ). 



Based on the framework of our previous modeling (lHan et al.ll2009i ). we have constructed new 
synthetic CMDs for NGC 1851. The best fitting parameters were determined with the technique 
described in §2.2. In the case of NGC 1851, the re is still no general consensus on the observed value 



of Af CNO /Fe] between th e two su bpopulations (|Yong et al.ll2009i : IVillanova et al.ll2010l : lLardo et al 



2O120 . but iGratton et al.l (|2012bl ) recently concluded that the difference in [CNO/Fe] is less than 
0.2 dex. Considering this and the similarity of NGC 1851 with M22, here, we have adopted 0.1 dex 
for A [CNO/Fe]. Given the small value, this should have only a little effect in our modeling. 

In Figure 14, our models are compared with new CMDs of NGC 1851, obtained from Ca-by 
photometry (S.-I. Han et al. 2012, in pr eparation) . This obser vation confirms the previously re- 
ported split of the RGB in the hk index (|J.-W. Lee et al.ll2009allb|). Sim ilarly, Figure 15 compares 



our models with the CMDs from the ACS survey ( Milone et al.l I2008I ). 5 The model parameters 
determined by our fitting technique are listed in Table 4. With the small variations in metallicity 
(A [Fe/H] = 0.13 dex) and the total CNO content (A[CNO/Fe] = 0.1 dex) between the two sub- 
populations, our fitting technique suggests that the difference in helium abundance is considerable 
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(AY = 0.06 ± 0.01), while the difference in age is relatively small (At = 0.5 ± 0.2 Gyr). This 
result is qualitatively similar to the case of M22. The faint SGB can be explained mostly by the 
combined effects of the metallicity and helium enhancements, as both of them make MSTOs fainter 
(Fig. 15), while the small variations in the total CNO abundance and age have only little impact 
on the CMD. In particular, as discussed above (Fig. 7), the total CNO content and age basically 
affect MSTOs in opposite directions. If there were no difference in [CNO/Fe] between Gl and G2, 
then the age difference would be even smaller (At ~ 0.0 Gyr). The Ca-strong redder RGB in the hk 
index can be accounted for solely by the effect of metal enhancement, which makes RGBs redder 
(bottom panels of Fig. 14). The helium enhancement slightly moves the RGB to the opposite 
direction, but the effect of metal enrichment is far more important. In the HB region, the effect of 
helium enhancement, which makes HB morphology bluer, easily overcomes the effects of metallicity 
and CNO contents (see Fig. 13), and therefore the presence of blue HB, populated by G2, can 
be explained. Note that our model for G2 is also c onsistent with recent measurement of helium 
abundances for the blue HB stars (Y = 0.29 ± 0.05; faratton et aj]|2012bl ). 



Our models indicate, therefore, that the majority subpopulation, Gl, comprises the metal- 
poor bright SGB, Ca-weak bluer RGB, and red HB stars; while the minority subpopulation, G2, 
is associated with the metal-rich faint SGB, Ca-strong redder RGB, and blue HB stars. This is 
consistent with the population ratios estimated from Figures 14 and 15. The number fraction 
between the bright and faint SGBs is roughly 0.68:0.32 (Fig. 15c), and the ratio between the bluer 
and redder RGBs is about 0.71:0.29 ( Fig. 1 4c). This also agrees with the number ratios reported 
in previous works (jMilone et alJhood . 12003 ; |j.-W. Lee et al.ll2009al : lHan et alJhood ). to within the 
uncertainty. The ratio between the red and blue HBs is about 0.69:0.31 (Fig. 14a and Fig. 15a), 
if RR Lyrae stars belong to metal-poor subpopulation (Gl). Note that the mean period of RR 
Lyrae stars (<P/> = 0.531 day), which is normal for the metallicity of NGC 1851 (IWalkerlll998l : 



Clement et al.ll200ll ). also indicates that most of the RR Lyrae variables come from the majority 



subpopulation (Gl) having normal helium abundance. 



4. DISCUSSION 

We have compared our models with the observations of three peculiar GCs, u Cen, M22, 
and NGC 1851, to estimate ages, metallicities, and helium abundances of their subpopulations. 
Based on our results, in Figure 16, we schematically illustrate the star formation histories and 
accompanying enrichments of metallicity and helium abundance in these stellar systems. We find 
that the age differences between the metal-rich and metal-poor subpopulations in these GCs are 
relatively small, ~0.3-1.7 Gyr, and metal-rich subpopulations with redder RGBs are also enhanced 
in helium abundance. The formation time scale of stellar populations in these GCs is therefore 
expected to be fairly short, i.e., ~1.7 Gyr for oj Cen and less than 1 Gyr for M22 and NGC 1851. 

How this could be understood in terms of their formation origin? The most natural assump- 
tion would be that they are self-enriched, in which later generations of stars formed from the gas 
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enriched by the ejecta of massive stars in earlier generations. As for the sources of chemical enrich- 
ment and pollution in these GCs, at least three viable mechanisms have bee n proposed thus far . 
The explosion of type II SNe is required for the enrich ments of heavy elements (jTimmes et al.l l 1995) 
and perhaps helium ()Norrisl 12004 ; iPiotto et al.ll2005l ). while winds from the FRM Ss and IMAGBs 



are suggested for the pollution of light elements and the en hancement of helium (jDecressin et al 



2007 



D'Antona fc Ventural 120071 : 1 Ventura & D'Antonal l2008). The enhancement of the total CNO 



abundance in the metal- rich later generation subpopulations would also indicate the contribution 



by type II SNe (see also iMarino et al.l l2012al ) . Considering the order of time required to emerge 
for these mechanisms, a possible scenario would be as follows: (1) Formation of the first-generation 
metal-poor (bluer RGB) stars from the gas having normal helium and light element abundances. 
(2) Pollution of the remaining gas by the winds from FRMSs, which enhance helium and alter 
the abundance profile of light elements. (3) The most massive (> 8 M ) stars then would ex- 
plode as type II SNe, enriching overall metallicity, including heavy elements and the total CNO 
content. (4) Further pollution of the gas by the ejecta from IMAGB stars (~3-7 M Q ) would fol- 
low, which enhance helium and change light element abundances. (5) Finally, formation of the 
second-generation metal-rich (redder RGB) stars from the gas now enriched in overall metallicity 
and helium and polluted in light element abundances. While this scenario can naturally explain 
the presence of metal-rich and helium-enhanced subpopulations in these GCs, this requires that 
they were much more massive in the past, because their present masses are to o small to retain the 



ejecta of SNe explosions (e.g.. iDopita fc Smithlll98q ; iBaumgardt et al.ll2008l ). This suggests that 



Galaxy, as is widely accepted for oj Cen (e.g 



2005; lMeza et al 
2009 



2005 



Pflamm-Altenburg Kroupal 2009; 



Lee et a" 



Lee et alJl2007l:l J.-W. Lee et al 



Cohen et al 



1999; 



2009a; 



2010 



Bekki & Freeman 


2003; 


Piotto et al. 


Baumgardt et al. 


2008: 


?erraro et al. 


Bekki 


2012, and references therein). 



Note that a significant fraction (~30%) of the helium-enhanced subpopulation in these GCs is also 
best explained if the second generation stars were forme d from enriched gas t rapped in the deep 
gravitational potential well of the ancient dwarf galaxies (jBekki fc Norrisl 120061 ) . 



Recently, ICarretta et all (|2010cl . l2011bh and Bekki Yond (|2012j) suggested that NGC 1851 
might have formed by simple merging of two GCs having different heavy element abundances 
initially belonged to a proto-dwarf galaxy. While this scenario can explain the presence of RGB 
split and the spread in heavy element abundance without assuming self-enrichment, this is not 
easily reconciled with one of our main results, which suggests that the metal-rich subpopulation is 
also enhanced in helium abundance. In this scenario, it would be difficult to understand why all 
stars in initially more metal-rich GCs were selectively enhanced in helium abundance, while those 
in metal-poor GCs were not. Note, however, that this merging scenario might be able to explain the 
absence of helium enhancement between Gl and G2 in uj Cen (see §3.1), where Gl and G2 formed 
by simple merging of two GCs in the central part of a proto-dwarf galaxy, while more metal-rich 
and helium-enhanced subpopulations, G3, G4, and G5, formed later from the enriched gas trapped 
in the core of this merged stellar system in a proto-dwarf galaxy. Despite many uncertainties, all 
of the above scenarios are possible only in the primeval dwarf galaxy environment, and therefore, 
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dwarf galaxy origin of these GCs appears to be inevitable. 

We are indebted to the anonymous referee for encouraging us to include the CNO enhanced 
models and to employ parameter estimation technique, which led to several key improvements 
in the manuscript. We also thank Sang-Il Han for providing the observational data in advance 
of publication, Chongsam Na for his assistance in the construction of nitrogen enhanced stellar 
evolutionary tracks, and Suk-Jin Yoon for helpful discussions. Support for this work was provided 
by the National Research Foundation of Korea to the Center for Galaxy Evolution Research. 
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Table 1. Model Ingredients 



Ingredient References 



Isochrones from MS to RGB Yonsei-Yale (Y 2 ) isochrones (Yi et al. 2008 ; Y.-W. Lee et al. 2012, in preparation) 

HB evolutionary tracks Yonsei-Yale (Y 2 ) HB tracks (Y.-W. Lee et al. 2012, in preparation) 

Color-T c ff transformation . . ^reen_e^_aL Q987) color table 

Castelli & Kurucz Q2003) model atmosphere 
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Table 2. Parameters suggested from our best simulation of uj Cen 



Population Z [Fc/H] a A[CNO/Fe] b Y Age(Gyr) Mass-loss(M ) c Fraction d other name d 



Gl 0.0005 -1.81 0.0 0.231 13.1 ± 0.2 0.199 0.42 RGB-MP 

G2 0.0009 -1.55 0.0 0.232 13.0 ± 0.3 0.217 0.28 RGB-MIntl 

G3 0.0015 -1.31 0.14 0.41 ± 0.02 12.0 ± 0.4 0.171 0.17 RGB-MInt2 

G4 0.0057 -1.01 0.47 0.38 ± 0.02 11.4 ± 0.4 0.220 0.08 RGB-MInt3 

G5 0.0136 -0.62 0.47 0.39 ± 0.02 11.4 ± 0.5 0.251 0.05 RGB-a 



a [a/Fe] =0.3 

b From lMarino et al.l ll2012ah . 

c Mean mass-loss on the RGB for r\ = 0.53. 

d From lSollima et al.l d2005ah and lPancino et al.l ll2000h . 
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Table 3. Parameters suggested from our best simulation of M22 

Population Z [Fe/H] a A[CNO/Fc] b Y Age(Gyr) Mass-loss(A/ Q ) c Fraction 

Gl 0.00035 -1.96 0.0 0.231 12.8 ± 0.2 0.189 0.7 

G2 0.00067 -1.71 0.13 0.32 ± 0.04 12.5 ± 0.4 0.192 0.3 



a [a/Fe] = 0.3 

b From lMarino et all (|2012bl ). 

c Mean mass-loss on the RGB for r\ 



= 0.53. 



Table 4. Parameters suggested from our best simulation of NGC 1851 



Population 


Z 


[Fe/H] a 


A[CNO/Fo] b 


Y 


Age(Gyr) 


Mass-loss(M Q ) c 


Fraction 


Gl 


0.0012 


-1.43 


0.0 


0.232 


11.1 ± 0.1 


0.208 


0.7 


G2 


0.0018 


-1.30 


0.10 


0.29 ± 0.01 


10.6 ± 0.2 


0.212 


0.3 



a [a/Fe] = 0.3 

From Villanova et al. |2010h and lGratton et alj (|2012bl ) 

c Mean mass- loss on the RGB for r\ — 0.53. 
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Fig. 1. — Calibration of the iReimera (|1977l ) mass-loss coefficient, p, on the HB morpho logy vs. 
[Fe/H] relation. Data for the inner halo (R GC < 8 kpc) GCs from Lee et all f)l994j . 120071 ). Solid 
lines are theoretical isochrones produced by our synthetic HB models for 13 Gyr using rj = 0.46, 
0.53, and 0.60, respectively. Two metal-rich GCs with extended HB are marked and are excluded 
in the comparison. 
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Fig. 2. — (a)-(d) Comparison of iso chrones and ZAHB models constructed with the color-T e g- 
transformations by iGreen et al.l (|1987l . solid lines) and by ICastelli fc Kuruca (|2003l . dotted lines) . 
The model lines plotted are identical to the models for the most metal-poor (purple lines) and the 
most metal-rich (red lines) subpopulations in u Cen (see §3.1 below), (e) Solid lines are models 
corrected by a dding the difference betw een the two color transformations in (d) to the dotted lines 
obtained from ICastelli fc Kuruca (120031 ) . 
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(a) NGC 1851 Obs. 




0.45 0.50 0.55 0.60 0.65 0.45 0.50 0.55 0.60 0.65 
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Fig. 3. — An example of the comparison between the observed and synthetic Hess diagrams used 
in our \ 2 minimization tech nique for NGC 185 1. (a) Red points represent CMD regions used in our 
X 2 minimization (data from lMilone et aDbood ). (b)-(e) The observed and synthetic Hess diagrams 
for HB and SGB regions, respectively. 
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Age [Gyr] Age [Gyr] 

Fig. 4. — Distributions of reduced x 2 {xt) values in the fitting of our models with the observed 
Hess diagram for NGC 1851. For the first-generation (Gl) population, only age is examined, while 
for the second-generation (G2) population, both age and helium abundance (Y) are simultaneously 
investigated to find the values with minimum % 2 . In panel (b), the values of x 2 are listed on the 
contour map (see text). 
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Fig. 5. — Comparison of our models with the observations for oj Cen. (a), (c) Observed 
CMDs by HST A CS/WFC for the F435W, F625W and F606W, F814W passbands (data from 
Bellini et al.1 1201CI ). (b), (d) Our population models compared on the observed CMDs. Param- 
eters suggested from our best simulation are listed in Table 4. Adopted distant modulus and 
reddening are (m-M) F62 5W = 14.03, E(F435W-F625W) = 0.225, and (m-M) Fmm = 14.05, 
E(F606W-F8UW) = 0.11, respectively. 
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Fig. 6.— (a), (b) Sa me as Fig. 5 but for the F275W (HST WFC3/UVIS) and F814W (ACS/WFC) 
passbands (data from lBellini et al.ll2010l ). We can see that the model parameters obtained from the 
optical (F435W, F625W) bands can also reproduce the observed features from the UV (F275W) 
and broad I (F814W) bands, (c) CMD for the MS and SGB region in the ACS/WFC F435W and 
F658N passbands. (d) Our synthetic models constructed with the photometric errors and total 
number of stars comparable to the observed CMD in panel (c). Adopted distant modulus and 
reddening are {m-M) F27 5w = 14.85, E(F275W-F8UW) = 0.91, and (m-M) F658A r = 13.97, 
E(F435W-F658N) = 0.22, respectively. 
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Fig. 7. — Similar to Fig. 5b, but here the models are constructed to illustrate the effects of 
A[CNO/Fe] on the age dating, (a) If A[CNO/Fe] = 0.0 between Gl and G5, no difference in age 
is obtained between them, (b), (c) When the observed A[CNO/Fe] is included, At = 1.7 Gyr is 
required between Gl and G5, in the sense that G5 is younger, otherwise models would not match 
the observations. 
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Fig. 8. — Similar to Fig. 5b, but here the models are constructed to illustrate, independently, the 
effects of metallicity, the total CNO content, age, and helium abundance on the overall morphology 
in the CMD (see text). 
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Fig. 9. — The effect of helium enhancement on the stellar spectra, (a) Synthetic spectra 
(jCastelli Kuruca 120031 ) of two stars near the MSTO, with the same T e g = 6000 K, log g = 
4.5, and [Fe/H] = —1.5, but different helium abundance (Y = 0.25 & 0.35). Green lines are filter 
throughputs for three UV-passbands (F225W, F275W, and F336W) of HST WFC 3/UVIS. (b) 



Obse rved CMD of u Cen in F225W and ACS/WFC F814W passbands (data from iBellini et al 



2010) . (c) The effect of helium-enhanced spectra on the CMD. Color codings are as in Fig. 5. 
The dotted lines are models from the normal helium (Y = 0.25) spectra, while solid lines are those 
from the helium-enhanced (Y = 0.38—0.41) spectra. Adopted distant modulus and reddening are 
(m-M) FS uw = 13.9 and E(F225W-F8UW) = 1.26, respectively. 
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Fig. 10. — Comparison of the observed RR Lyrae variables (squares) in to Cen with the synthetic 
HB models (c r osses) . Re d squares are metal-ri ch RRab stars ([Fe/H] > —1.4). RR Lyrae data from 
Sollima et al.l (|2006l ) and iKaluznv et al.l (|2004l ) . Purple and blue crosses are model RR Lyrae stars 
from Gl and G2, respectively. This suggests that there is an additional minor subpopulation be- 
tween G2 and G3, having normal helium-abundance and [Fe/H] ~ —1.3. Adopted distant modulus 
and reddening are {m—M)y = 14.1 and E(B—V) = 0.14, respectively (see text). 
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Fig. 11. — Comparison of our models with the observations for M22. (a), (c) Stromgren b, y, 
and calcium {Co) narrow-band photometry from S.-I. Han et al. (2012, in preparation), where the 
hk index is defined as hk = (Ca — b) — (b — y). (b), (d) Our population models compared on the 
observed CMDs. Adopted distant modulus and reddening are (m—M) y = 13.80, E(b—y) = 0.255, 
and E(hk) = -0.012. 
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Fig. 12. — Comparison of our models with the observations for M22. (a) Observed CMD by 
HST ACS/WFC for the F606W and F814W passbands. Photometric data from the ACS survey 
(jSarajedini et al.1 120071 : lAnderson et al.ll2008l ). (b) Our models compared on the observed CMD. 
(c), (d) Same as (a) k, (b) but zo omed in for the SGB region. A clearer view of the SGB split 
can be seen in lMilone et al.l (|2010l . see their Fig. 2). Adopted distant modulus and reddening are 
(m-M) F60 Q W = 13.69 and E(F606W-F8UW) = 0.355, respectively. 
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Fig. 13. — The effects of variations in the total CNO abundance, metallicity, and helium abundance 
in our models for M22 (see text). 




Fig. 14. — Same as Fig. 11, but for NGC 1851 (data from S.-I. Han et al. 2012, in preparation). 
Adopted distant modulus and reddening are {m—M) y = 15.62, E(b—y) = 0.015, and E{hk) = 
-0.015. 
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Fig. 15.— Same as Fig. 12, but for NGC 1 851. Data from the A CS survey (jMilone et alJl2008l ). A 
clearer view of the SGB split can be seen in lMilone et al.l (J201CJ, see their Fig. 2). Adopted distant 
modulus and reddening are (m— M)fqqqw = 15.52 and E(F60QW— F814W) = 0.025, respectively. 
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Fig. 16. — Schematic illustrations for the star formation histories and accompanying enrichments 
of metallicity and helium abundance in three GCs considered in this paper. The formation time 
scales of all subpopulations appear to be ~1.7 Gyr for uj Cen and less than 1 Gyr for M22 and 
NGC 1851. Metal-rich subpopulations in these GCs are also enhanced in helium abundance. There 
is a dramatic increase of helium abundance between G2 and G3 in uj Cen. One sigma error bars 
for age and helium abundance are from our x 2 minimization, while the uncertainty in [Fe/H] is 
assumed to be 0.1 dex. 



